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SOLITARY GRAVITY-CAPILLARY WATER WAVES WITH 

POINT VORTICES 


KRISTOFFER VARHOLM 


Abstract. We construct small-amplitude solitary traveling gravity¬ 
capillary water waves with a finite number of point vortices along a 
vertical line, on finite depth. This is done using a local bifurcation 
argument. The properties of the resulting waves are also examined: We 
find that they depend significantly on the position of the point vortices 
in the water column. 


1. Introduction 

The steady water-wave problem concerns two-dimensional water waves 
propagating with constant velocity and without change of shape. Historically, 
the focus has mainly been on irrotational waves, which are waves where the 
vorticity 1 

uj := V x w = v x — u y 

of the velocity held w = (it, v) is identically zero. One reason for this is 
Kelvin’s circulation theorem [29,31], which says that a flow which is initially 
irrotational will remain so for all time, as long as it is only affected by 
conservative body forces (e.g. gravity). Another reason is mathematical, 
as the velocity held can then be written as the gradient of a harmonic 
function; the velocity potential. This enables the use of powerful tools from 
complex- and harmonic analysis, and the problem can be reduced to one on 
the boundary in a number of different ways [1,38]. An important class of 
such waves are the Stokes waves, which are periodic waves that rise and fall 
exactly once every minimal period. The Stokes conjecture on the nature of the 
so-called Stokes wave of greatest height fueled research on waves throughout 
the 20th century, and would not be fully resolved until 2004 (see the survey 
[43] and [39], which settled the convexity of this wave). 

More recently, however, there has been renewed interest in rotational 
waves. There are several situations where such waves are appropriate, as 
effects like wind, temperature or salinity gradients can all induce rotation 
[37]. Rotational waves can be markedly different from irrotational waves: For 
instance, in rotational waves it is possible to have internal stagnation points 
and critical layers of closed streamlines known as cat’s eye vortices [15]. 
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The first result on rotational waves came surprisingly early, in the beginning 
of the 1800s with [22] (for a more modern exposition, see [4]). There, Gerstner 
gave the first, and still the only known, explicit (nontrivial) gravity-wave 
solution to the Euler equations on infinite depth. Although significant because 
it is an exact solution, it is viewed as more of a mathematical curiosity, even 
today (see [5, Chapter 4.3]). Much later, in [14], came the first existence 
result for small-amplitude waves with quite general vorticity distributions. A 
vorticity distribution is a function 7 : R —> R such that 

= 7 (Y0> 

where i/j is the relative stream function (which, unlike the velocity potential, 
is still available for rotational waves, but is not harmonic). A sufficient, but 
not necessary, condition for such a vorticity distribution to exist is that the 
wave has no stagnation points. Several improvements have been made to the 
existence result of Dubreil-Jacotin, but it was not until the pioneering article 
[7] that large waves were constructed, using global bifurcation theory. This 
article sparked mathematical research into rotational waves. 

The use of a semi-hodograph transform in [14] and [7], and the correspond¬ 
ing deep-water result in [27], means that the resulting waves cannot exhibit 
critical layers. Since then, small-amplitude waves with constant vorticity 
and a critical layer have been constructed in [46], and later in [9] with a 
different approach that allows for waves with overhanging surface profiles 
(there is numerical evidence for the existence of such waves, e.g. [44], but this 
is still an open problem). A reasonable next step is that of waves with an 
affine vorticity distribution, whose existence was shown in [16,19]. Spurred 
by the above results there has also been interest in studying the properties 
and dynamics of these waves below the surface [15,46]. This had been done 
for linear waves in [18]. Several other avenues have also been considered: We 
mention heterogeneous waves both with [25,48] and without [20,47] surface 
tension, waves with discontinuous vorticity [ 8 ], a variational approach [ 2 ] and 
Hamiltonian formulation with center manifold reduction [24]. Existence of 
large amplitude waves with constant vorticity and a critical layer was estab¬ 
lished in [36], in the presence of capillary effects. There is also a forthcoming 
result for pure gravity waves [3], using an entirely different approach. 

Common for all the previously mentioned works on rotational waves is 
the feature that the vorticity is supported on the entire fluid domain (due 
to the assumption of the existence of a vorticity distribution). Recently, 
gravity-capillary waves with compactly supported vorticity were constructed 
in [41], on infinite depth. This includes small- and large-amplitude periodic 
waves with a point vortex, and small-amplitude solitary waves with either a 
point vortex or vortex patch. By a point vortex we mean that the vorticity is 
given by a (5-function, while we use vortex patch to mean that the vorticity is 
locally integrable and compactly supported. The waves with a point vortex 
are the simplest form of waves with compactly supported vorticity, and are 
in a sense “almost irrotational”. 

In this paper, which is based on [45], we extend the existence result for 
solitary small-amplitude waves with a point vortex to finite depth, and also 
give both qualitative and quantitative properties for these waves. The main 
approach to showing existence follows that of [41], but we also treat the 
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natural generalization of waves with several point vortices along a vertical 
line, and show existence for all but exceptional configurations of vortices. 
Finally, by finding an explicit expression for the rotational part of the stream 
function, we give some explicit expressions for the small-amplitude periodic 
waves with a point vortex on infinite depth that were constructed in [41]. 

An outline of the paper is as follows: In Section 2 we formulate the problem, 
and in Section 3 we give the functional-analytic setting for this formulation. 
Then, in Section 4 we prove existence of small solutions, and give some 
properties for these. Section 5 treats the extension to several point vortices. 
The final section, Section 6, contains the explicit expressions for periodic 
waves on infinite depth. 


2. Formulation 

Under the assumption of inviscid (absence of viscosity) and incompressible 
(constant fluid density) flow, the governing equations of motion are the so- 
called incompressible Euler equations. For describing water waves on the 
open sea, these are realistic assumptions [29,31], and standard. We will 
further assume two-dimensional flow under the influence of gravity, where 
the Cartesian coordinates (x, y ) describe the horizontal and vertical direction, 
respectively. Then the equations read 

u>t + (w ■ X7)w = —Vp — ge 2 , (Conservation of momentum) 

V • w = 0, (Conservation of mass) 

where w = (u, v ) is the velocity of the fluid, p is the pressure distribution 
and —ge 2 = (0, —g) is the constant gravitational acceleration 2 . 

For convenience we place, at time t. the flat bottom at 

{(x,y) e K 2 : y = -h} 

and the surface at 

{(x,y) G M 2 : y = r)(x,t)}, 

where 77 describes the deviation of the free boundary. We assume that 
is bounded, continuous and strictly bounded below by —h. It should be 
emphasized that, due to the free boundary assumption, the function p is a 
priori unknown; determining it is part of the problem. 

In addition to Equation (2.1), we require boundary conditions to match our 
domain. In order to model the bottom being impermeable, we will demand 
that 


v\ y= _ h = 0 (Kinematic boundary condition at bottom), 

with which we mean that v{x, — h,t ) = 0 for all x and t. Next, we impose 
the condition 

ur) x + rjt = v (Kinematic boundary condition at surface) 


2 The constant g is approximately 9.8m/s 2 , varying by less than 0.4% on the Earth’s 
suface (see [26]). 
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at the surface. This equation is what connects the free boundary to the fluid, 
and is equivalent to demanding that particles at the surface will remain there. 
We also require that 

(2.2) P\y=rj = a 2 K(?7), (Dynamic boundary condition) 

where a 2 > 0 describes the surface tension and k is the nonlinear differential 
operator 

Vx \ Vxx 


k(v) ■= 


\3’ 


MJ (vx 

yielding the curvature of the surface. The symbol (• 


denotes the Japanese 


bracket defined through x H > (1 + \x I 2 ) 1 / 2 . Equation (2.2) is known to 
physicists as the Young-Laplace equation, and states that the pressure 
difference across a fluid interface (in this case water/air) is proportional to 
its curvature. 

Note that in the lower limit, a 2 = 0, the dynamic boundary condition 
in Equation (2.2) corresponds to the assumption of constant pressure on 
the surface, but we will require that a 2 be strictly positive. The proof of 
Theorem 4.1, for example, relies upon the assumption that a 2 > 0. 


2.1. The vorticity equation. By taking the curl of Equation (2.1), one 
obtains after some simple calculations that 

(2.3) ut + V • (cjw) = 0, 

which states that the vorticity oj is transported by the vector field w. Due 
to this, it is natural to expect that if the vorticity consists of a point vortex 
at some time, then it will remain a point vortex at all future times, and be 
transported with the flow. It should be emphasized that, for now, this is 
not justified by Equation (2.3); the multiplication of u with w is not well 
defined, as w will not be smooth at the point vortex. Thus, we will have need 
of a weaker form of the equation. We remind the reader of the fundamental 
solution of the Poisson equation. 

Proposition 2.1 (Newtonian potential). The distribution T 6 Lj( ic (M 2 ) 
defined by 

T(x,y) := ^log {x 2 + y 2 ) 

satisfies 

v 1 r ( i, ! /) := 

and 

AT = V x V^T = 6. 


If uj is of the form 

Ld(t) = 8(xo(t),y 0 (t))i 

then we deduce from Proposition 2.1 that w is of the form 


w(x, y, t) 


i (yo(t) - y,x - x 0 (t)) 

2 vr (x - x 0 (t)) 2 + {y- yo{t)) 2 


+ w(x,y,t), 


where w satisfies V • w = 0 and V X w = 0, and is therefore smooth in space 
(see the discussion before Equation (2.10)). As the first term, which we may 
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think of as the part of w generated by the point vortex, is singular, divergence 
free and odd around (xo(t),yo(t)), it is not unreasonable to think that the 
dynamics of the point vortex should depend only on w. In other words, that 
the path t H > (xo(t),yo(t)) along which the point vortex moves should satisfy 

(2.4) (±o, yo) = w. 

This can indeed be made rigorous. In [32, Theorems 4.1 and 4.2] it is proved 
that if one considers initial data consisting of a vortex patch converging in the 
sense of distributions to a point vortex, then the weak solutions of the vorticity 
transport equation converge to a moving point vortex in an appropriate sense. 
Moreover, the position of this point vortex satisfies Equation (2.4). Thus, we 
will allow for point vortices, as long as they are propagated in the fluid as in 
Equation (2.4). 


2.2. Traveling waves. We now assume that there are functions w,p,fj, 
depending only on space, and a constant velocity c £ M such that 

w(x, y, t) = w(x - ct,y), 

P(x,y,t) = p(x-ct,y), 
p(x, t ) = fj(x — ct ) 

for all relevant x. y and t. Positive and negative c then correspond to waves 
moving in the positive and negative x-directions, respectively. In the new 
steady variables (x,y) = (x — ct,y ), after dropping the tildes, our equations 
read 


(Conservation of momentum) 
(Conservation of mass) 


at y = —h, 
at y = T](x), 


(Kinematic) 

(Kinematic) 

(Dynamic) 


(2.5) (w ■ S/)w — cw x = —'S/p — ge 2 , 

V • w = 0, 

with boundary conditions 

(2.6) v = 0, 

(2.7) (■ u — c)rj' = v, 

(2.8) p = —a 2 n(p), 

on the now time-independent domain 

D(?y) := {(x,y) £ M 2 : —h < x < g(x)}. 

We call the problem of finding w , p and p such that these equations are 
satisfied the steady water-wave problem. Note also that the vorticity equation 
given in Equation (2.4) reduces to 

(2.9) (c,0) = w(x 0 ,yo) 
for a point vortex centered at (xq, ]Jq) £ f2(r/). 


2.3. The Zakharov—Craig—Sulem formulation. It turns out that it is 
possible to reduce the water-wave problem to an entirely one-dimensional one 
on the surface in a clever way. This is known as the Zakharov-Craig-Sulem 
formulation, and was first introduced by Zakharov in [49], and then later put 
on a firmer mathematical basis in [10,11]. The original formulation relies on 
the fluid being irrotational, but it is in fact sufficient that this holds near the 
surface. This is where the compact support of the vorticity comes in. 



6 


KRISTOFFER VARHOLM 


Suppose that we have solved the steady water wave problem for some 
w,p,rj. It is then convenient to split the velocity w as 

w = w + W, 

where w is irrotational, that is, V X vj = 0 and V x W = uj. We also assume 
that both w and W are divergence free. For us, W will be known. 

Although we will allow for uj to be a singular distribution, the vector field 
w = (ft, v) will be assumed to at least be in the Sobolev space i7 1 (fl(r/)) 2 , and 
W = ( U , V) to be at least L\ oc ({\(r])). By the assumption of \7-u> = V-W = 0, 
the differentials 

v dx — u dy , V dx — U dy 

on £l(r)) are closed. Hence, as 12 (r/) is simply connected, these differentials 
are exact by generalizations of the Poincare lemma [33, Theorems 2.1 and 
3.1]. Thus, there are stream functions t/> £ , T £ Lj oc (Q(r])) (by 

[13, Corollary 2.1]), determined uniquely modulo constants by w and W, 
such that 

w = V ± 'ip, W = V ± T. 

Moreover, by the assumed curl of these vector fields, the function ^ is 
harmonic, and T satisfies AT = uj. In particular, i/j is smooth, and so is W 
outside the support of uj. 

By the above, we thus have that 

(2.10) w = V " 1 {$> + T) 

holds for the velocity held w. Suppose now that W and T are chosen such 
that T = 0 at the bottom. Then the boundary condition in Equation (2.6) 
translates to ift being constant along the bottom. Since is unique modulo 
constants, we may as well take this condition to be 

=0 

y=—h 

instead. 

We will now apply the assumption that suppw C 12 ( 77 ) is compact. This 
means that a velocity potential for w exists on any simply connected domain 
not containing supp uj. Equation (2.5) can then in turn be used to show that 

V cu+ 77 M 2 +P + gy^j = 0 

holds on any such domain. Hence, in particular, we obtain the surface 
Bernoulli equation 

( 2 . 11 ) c(4>y + Ty) + ^|Vt/> + VT | 2 - a 2 n(r]) + gy = C, at y = rj(x), 

where C is a real constant. Here we have inserted the boundary condition 
for the pressure at the surface, Equation (2.8). We will set (7 = 0, since 
we are looking for localized waves. This can be seen by letting |x| —> 00 in 
Equation (2.11). We now need the following formal definitions to proceed, 
which will be specified later on. 
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Definition 2.2 (Harmonic extension operator). Given 77, we define the 
harmonic extension operator H (77) as the operator mapping a function £: M —> 
M to the harmonic function ip: H(t/) —> M satisfying 

= C> 

= 0 . 

Definition 2.3 (Dirichlet-to-Neumann operator). Given r 7, we define the 
Dirichlet-to-Neumann operator G(jj) as the operator mapping Dirichlet data 
to non-normalized Neumann data; that is, the operator defined by 

G(t7)C:= (- V ',l)-V[H( V K\\y= v 

for functions (: M —> M. 

With Definition 2.2 in mind, define (: M -> M by 

(2-12) C 

that is, the trace of ip on the surface. By our assumptions, then, we have 

$ = H(r])C, 

and we will use this to reformulate Equation (2.11) in a way that only involves 
£ and ’E. Note that 

C' = i’x+v'i’y, 

G(r ?)C = ~r\^ x + ip y , 

where the right-hand side is evaluated at y = r]{x). Inverting these relations 
and inserting them in Equation (2.11) yields 

KC' + g(t?)c 

L W) 

(2.13) (C' + (1, rf) • V^) 2 + (G(t 7)C + (-r/, 1) • Wh) 2 

2(7/> 2 

+ gij — a 2 n(r]) = 0, 

and in a similar fashion, one obtains 

(2.14) c?/ + C' + (I,?/) • V'h = 0 

from the kinematic boundary condition in Equation (2.7). Equation (2.14) 
can be integrated once to yield 

(2.15) cr]+ ( + ^=0, 

where we have assumed decay at infinity. We emphasize that the function 4' 
and its derivatives are evaluated at y = ij(x) in Equations (2.13) to (2.15), 
which we suppress for readability. Equations (2.13) and (2.15) form the 
Zakharov-Craig-Sulem formulation, which we will combine with a suitable 
vorticity equation. One may note that the pressure, p, has been eliminated 
from the formulation entirely. 

Remark 2.4. Equation (2.13) is slightly different than the equation used 
in [41]. The equation in [41] can be obtained by inserting the kinematic 
boundary condition, Equation (2.14), into Equation (2.13). 
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3. Functional-analytic setting 

We now focus on proving the existence of a family of small amplitude 
and small velocity traveling waves with vorticity consisting of a point vortex 
situated on the y-axis. In other words, solutions with vorticity of the form 

u = e5 g , 


where 0 < |e| <C 1, 0 G (0,1) and where we have defined 

fie '■= fi(0,-(l-6)h)- 

The constant 9 then corresponds to the relative position of the point vortex 
above the bottom, and the parameter e, describing the strength of the vortex, 
will be used as the bifurcation parameter. 

We will from here on always assume that r/ is such that ?/(0) > —(1 — 9)h , 
which prevents the surface from touching the point vortex. Furthermore, we 
will also assume that y(0) < (1 — 9)h. The reason for this is purely technical 
(as we will see after Proposition 3.1). For the purpose of accounting for these 
assumptions, define the set 


(3.1) A e := { V G BC(R) : y > -h, | ?? ( 0 )| < (1 - 0)h}, 


whose intersection n H S (R) is open in H S (R) for any s > 1/2, by the 
Sobolev embedding H S (R) ^A BC(R). 

The next proposition describes the stream function that we will use for 
the rotational part of the velocity; the counterpart of the function G in [41]. 
While we could use a similar stream function on finite depth, it is more 
beneficial to work with one that is tailored for finite depth. 


Proposition 3.1 (Stream function). Let r] G Ag and define 4>: Q(r/) -A M by 


X >/ T v ) . = J_ w ( cosh(™/h) + cos(ir(y/h ~ 6)) \ 
47t \cosh.(irx/h) + cos(tt (y/h + 6)) J 

Then $ defines a regular distribution, and 


(3.2) 


Ad- = 5 g , 


$1 


y=0 




y=-h 


= 0 , 

= 0 . 


Moreover, the function (x, y) i-a 4>(x, y) — T(x, y + (1 — 9)h) is harmonic and 
satisfies 


(3.3) V x ($ - T(s, y + (1 - 0)h)) (0, -(1 - 6)h) 


4 h 


COt(7T0), 0 


where T is the Newtonian potential introduced in Proposition 2.1. 


Proof. We will apply Theorem A.l (see Appendix A) to prove this result, 
and thus need a bijective conformal map from the strip M x (—h, 0) C C to 
the unit disk D, mapping the point —i{ 1 — 9)h to the origin. This is done in 
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three steps: 


x (—h, 0) 

( z+ih)/h 

l x (0,1) 



exp(7T2:) 


x (0, oo) 


The conformal map for each individual step is well known from elementary 
complex analysis, see for instance [21, II.7 and p. 60]. Hence 

gK(z-\-ih)/h _ ginO 

f(z') * gir(z-\-ih)/h _ g—inO 


defines the desired map from the strip to the unit disk. By the aforementioned 
theorem, then, 


$(x,y) ■= — log(|/(x + iy)|) 


= iT og 


cosh(7rx//i) + cos(7r (y/h — 9)) 
cosh(7 rx/h) + cos(7r(y//i + 9 )) 


solves Equation (3.2) in M X (—h, 0). Because / extends to a meromorphic 
function on C, it is immediate that also solves Equation (3.2) in 11 (r/) 
(recall that ij £ Ag). 

Finally, we have 


V x ($-T(x,y + (l-(?)/»)) (0,(9) 


j (f"m\ 

4vr V f(i9) J 

(TcotM),o) 


by the final part of the same theorem, which will be important for the 
asymptotic velocity of the traveling waves that we shall obtain in Theorem 4.1. 

□ 


Note that the stream function 4> introduces a “mirror vortex” at (0, (1 — 9)h), 
and moreover is 2/r-periodic in the y-direction. This is the reason for the 
limitation on the height of the surface profiles in the set Ag defined in 
Equation (3.1). 

The next proposition is crucial, because the traces of $ and its derivatives 
on the surface enter in the Zakharov-Sulem-Craig formulation of the problem. 
Having an explicit expression for <J> enables us to prove the proposition in a 
quite direct way. 

Proposition 3.2. Suppose that rj £ H s (M) n Ag, where s > Then 

4>(-,y(-)) £i7 s (M), 

V ± 4>(-, ??(•)) € H S (R) 2 . 

Moreover, the dependence on rj is analytic. 
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Proof. We will only treat <F, as the argument for the derivative is similar. 
Observe that it is sufficient to consider the function defined by 

(3.4) x eA log(l + cos (rj(x) — 9) sech(x)) 

for r] £ H S (M) such that |i?(0)| < ir — 9, where 9 £ ( 0 , 7 r). Since 

cos(r/(x) — 9) sech(x) = (cos (r](x) — 9) — cos (9)) sech(.x) + 008 ( 6 *) sech(x) 

it follows by [40, Theorem 4 of 5.5.3], sech £ H S (R) and H S (R) being an 
algebra that the function in Equation (3.4) lies in H S (R ) and that the 
dependence on r] is analytic. Another application of the result in [40] then 
yields the desired result. □ 

As we have seen, because of the reliance on the stream function and the 
operators H(ij) and G(r]), a central problem is the solution of the Laplace 
equation, 

Af) = 0 in fl(r/), 

= C, =0 

y=v y=-h 

on the fluid domain 14(r/), given i] and £. We have the following theorem, which 
is adapted from Corollary 2.44 in [30], and which establishes both existence 
and uniqueness to Equation (3.5) in suitable Sobolev spaces. Functions on the 
surface will be identified with functions on the real line as in Equation (2.12). 

Theorem 3.3 (Well-posedness of the Laplace equation [30]). Suppose that 
r] £ H S (R )nAe for some s > 3/2, and that ( £ H 3//2 (R). Then Equation (3.5) 
has a unique solution in H 2 (14(?/)). 

Remark 3.4. While the natural setting for the velocity potential or the stream 
function on infinite depth is that of homogeneous Sobolev spaces, used in 
both [41] and [30], this is not the case for the stream function on finite 
depth. Because we require to be equal to a constant at the bottom, it 
must necessarily be the case that t/j tends to the same constant at infinity. 
Otherwise, because of the finite depth, tf y would not decay at infinity (in the 
sense that limu^, t y)\^oo ^yipi y) = 0 ), and therefore not describe a localized 
wave 3 . 

Theorem 3.3 enables us to define the harmonic extension operator described 
in Definition 2.2 as an operator H 3 ^ 2 (R) -A H 2 (Q(r])), and using this, defining 
the Dirichlet-to-Neumann operator. We refer the reader to [30], which is a 
rich source of results for these operators also in a more general setting. The 
results there are proved for the velocity potential, but should be adaptable 
for the stream function . 4 The below theorem is an amalgamation of parts 
from Corollary 2.40 and Theorems 3.15 and A.11 in [30]. See also [42], 

3 One could say that such a wave is localized if the limit exists and is different from zero, 
but this does not yield any new waves (only a change in the frame of reference). 

4 Compare with Theorems 3.49 and A. 13 in [30] for the case of infinite depth, where 
the boundary conditions for the Laplace equation for the stream function and velocity 
potential coincide. 
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Theorem 3.5 (Boundary operators [30]). Let s > 3/2 and suppose that 
r] G H S (R) n Ag. Then the harm.onic extension operator H(r]) and the 
Dirichlet-to-Neumann operator G(rj) are members of £?(.ff 3 / 2 (R), H 2 (Q(r])) 
and B(H S (R), respectively. The norms of these operators are 

uniformly bounded on subsets of H s (M) n A g that are bounded in the norm 
on H S (R). Moreover, the map G(-)( is analytic for fixed f G H S (R). 

In the same setting as in Theorem 3.5, the curvature of the surface is well 
defined: 

Proposition 3.6 (Curvature). The curvature operator k is well defined as an 
operator H S (R) —> H S ~ 2 (R) for any s > 3/2. Moreover, the map is analytic. 

Proof. Observe that the function /: M — > M defined by f(x) = x(x )~ 1 is 
smooth and satisfies /(0) = 0. As s — 1 > }, the result [40, Theorem 4 of 
5.5.3] ensures that f(r}') G PP _1 (M). Since / is also analytic, k is analytic by 
the same result. □ 

There is one thing we have not yet looked at, namely the vorticity equation 
Equation (2.9). Recalling Equation (2.10), we will consider velocity fields of 
the form 

(3.6) w = V x (H(rjX + e$) • 

From Proposition 2.1 we know that the part of the stream function that is 
generated by the point vortex at (0, —(1 — 6)h) is given by the Newtonian 
potential 

sT(x,y + (l - 0)h), 
whence Equation (3.6) reduces to 

(c,0) = V- L [iT(77)C](0, —(1 - 6)h) + e cot(7T0), 0^ , 

by Equation (3.3) in Proposition 3.1. 

In particular, this means that any solution necessarily must satisfy 

[H(vK]x(0,-{l~9)h) = 0. 

For simplicity, we choose to look for //, f in appropriately chosen subspaces of 
iP s (R), such that this condition is automatically satisfied. Specifically, define 

HLeniM) ■= {/ e H s (R) : / is even}, 

which is closed in H S (R), and therefore a Hilbert space in the inherited 
norm. We mention that it is still an open question whether asymmetric 
traveling waves exist. However, it is known that, in many situations, certain 
properties imply symmetry (see for instance [6,28]). Furthermore, under 
suitable assumptions, all symmetric waves are traveling waves [17]. 

Assume now that r] G Hf ven (R) n Ag, with s > 3/2, and that f G H^le. d(M). 
Then it must necessarily be the case that is even in x. Hence [H(r])£] x 

vanishes along the y -axis and so the vorticity equation reduces further to 

(3.7) c = c\£ — [H(r])f]y(0, — (1 — 0)h), where ci :=—cot(7r6 l ). 
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Observe also that the Dirichlet-to-Neumann operator G(rj) is well dehned as 
an operator H^ ven {R) -A for 77 G i/| ven n A q and s > 3/2, and that 

k can be viewed as an operator i7<f ven (M) —>• H^~ e „(R). 

Remark 3.7. One has to be careful with claims about the solution set when 
e = 0. Equation (3.7) of course actually only needs to be satisfied if e 7 ^ 0. 
This means that if we impose Equation (3.7), then we lose the trivial set of 
solutions 

(77,Cc,e) € { 0 } x { 0 } xRx { 0 } 

for the other equations, except for the point (0,0,0,0). This should be kept 
in mind in any claims of uniqueness. 

For convenience, define now the spaces 

X s := H! ven (R) X H s even (R) X R, 

Y* := X tflvenW X R 


and the set 


U^-.= {(r 1 ,C,c)eX s : V eA e }CX s , 


which accounts for the limitations on 77 . 

We proceed to introduce three maps that together will form the basis for 
our argument. For s > 3/2 we define F \: Ug X R —> by 


Fi{v,C,c,s) = c 


r/C_ + G(r])( 

W ) 2 


+ Ed 1 ), 


4- (C + ^(l.'/)W^ + ( G ( ? K + e(-,M).V1>f + _ j 


the map F2 : Ug X 


Hi 


m 2 

R) by 


F 2 {ri,(,c,e) = c ?7 + C + £-d>, 

and finally the map F 3 : Ug x R —> R by 

(3.8) F 3 (r],C,c ,£) = c- ci£+ [i7(ry)C] y (0, -(1 - 6 »)h). 

In all of these definitions, we really mean the traces of <f> and its derivatives 
on the surface. The pointwise evaluation in the second term of Equation (3.8) 
is allowed because is harmonic. It should be clear that all three maps 

F\. F 2 . F 3 are smooth. 

We can now define F : t/| X M —» Y s by 

F := (Fi,F 2 , F 3 ), 

and our task will then be to find solutions of the equation 


(3.9) F(r)X,c,e) = 0. 

One may immediately note that F( 0, 0,0, 0) = 0, so that the origin is a trivial 
solution. It will turn out that in a small neighborhood of the origin in X s X M, 
there is a unique curve of nontrivial solutions parametrized by the vortex 
strength parameter e. 
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4. Local bifurcation 

We can now finally state and prove the following theorem, establishing 
the existence of small, localized, traveling wave solutions with a point vortex. 
For this, we will use an implicit function theorem argument on F. Note that 
while we do not apply Crandall-Rabinowitz theorem [12], the situation is 
very much in the spirit of that theorem. We bifurcate from the family of 
trivial waves described in the remark after Equation (3.7) by introducing the 
vorticity equation. 

Theorem 4.1 (Traveling waves with a point vortex). Let s > 3/2 and 
0 £ (0,1). Then there exists an open interval I 3 0 and a C°°-curve 

I -> (H s even (R) n A 0 ) x H s even ( R) xRxl 

£ i y (r/(£),C(e),c(£),£) 

of solutions to the Zakharov-Craig-Sulem formulation, Equation (3.9), for a 
point vortex of strength e situated at (0, — (1 — 9)h). The solutions fulfil 

ry(£) = g 2 e 2 + 0(£ 4 ), 

(4.1) COO = C3& 3 + 0(e 4 ), 

c(e) = C\£ + C 3 £ 3 + 0(£ 4 ), 

in their respective spaces as e —> 0, where r/ 2 € i7f„ en (M) is defined by 
m ■= ~{g ~ a 2 dl)~ x x, X := ci$„(-, 0) + ^$„(-, °) 2 , 

and where 

C3 := i?2(ci + 4>y(-,0)), 

c 3 :=-[H(0)( 3 } y (0,-(l-e)h), 

with c\ as in Equation (3.7), <& as in Proposition 3.1 and H as in Defini¬ 
tion 2.2. 

Moreover, there is a neighborhood of the origin in t/| X M such that this 
curve describes all solutions to F(rjy,c,£) = 0 in that neighborhood. 

Proof. As remarked at the end of Section 3, the origin is a trivial solution. 
In order to apply the implicit function theorem, we require the first partial 
derivatives of F at this point. A direct calculation yields 

g — a 2 d 2 0 0 

(4.2) D X F( 0,0,0,0)= 0 0 > 

0 [77(0)4(0, -(l-6)h) l 

where the subscript X denotes the partial derivative with respect to the 
variable ( 77 , C, c) in X s . 

Now, every operator on the diagonal of D\F( 0, 0, 0,0) is an isomorphism. 
Indeed, the operator 

[g - a 2 d 2 x ] : 77 e s ven W ^ H s e ~ 2 n (R) 

corresponds to the Fourier multiplier g+d 2 ^ 2 . Since g , a 2 > 0, this operator is 
invertible, with inverse corresponding to the multiplier (i^+Q! 2 ^ 2 )^ 1 . The other 
two operators on the diagonal are identity operators, and therefore trivially 
invertible. Hence D\F( 0,0,0,0) £ B(X S ,Y S ) is also an isomorphism. 
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Thus we can use the implicit function theorem to conclude that there is 
an open interval I containing zero, an open set VC Uq containing ( 0 , 0 , 0 ), 
and a map / £ C°°(I , V) such that for ( 77 , £, c, e) £ V x /, we have 

F(v,C,c,s) = 0 «=>- (r?,C,c) = /(e). 


Furthermore, we obtain 


Df(0) = D x F[0, 0, ci, 0) — 1 Z) £ V(0, 0, ci, 0) 


0 

0 

ci 


which yields the first order terms in Equation (4.1). The higher-order terms 
can be obtained by inserting expansions for r](e), £(e) and c(e) into the 
equation F(rj(£), £(e), c(e), e) = 0. This concludes the proof of the theorem. 

□ 


Remark 4.2. Because Theorem 4.1 holds for any s > 3/2, we can get arbitrarily 
high regularity on the solutions, by possibly making the interval I smaller. 
We have not been able to conclude that they are smooth, however, since the 
interval could possibly shrink to a point as s —> oo. 


Observe that, because c i changes sign at 6 = 1/2, the direction in which 
the waves obtained in Theorem 4.1 will travel (for small e) depends on where 
the point vortex is in relation to the line y = —h/2. This does not come into 
play for waves on infinite depth. Note that if 0 = 1 — 1/h (when h > 1) then 

Cl = -^+°(l/h 2 ) 

as h —> oo, which is in agreement with what was found in [41] for a point 
vortex situated at ( 0 , — 1 ) on infinite depth. 

Since c\ vanishes when 9 = 1 / 2 , also the next term in the expansion for 
c(e) is of interest. We gave an expression for C 3 in Theorem 4.1, but have 
not determined its sign yet. We will treat the sign of C 3 after Theorem 4.6, 
which establishes some properties of the function re¬ 
written out, we have 


(4.3) 


X{x) = 


1 1 + COS(7T0) COsh(7 Tx/h) 


8 h 2 (cosh( 7 xx/h) + cos( 7 r 0)) 2 ’ 

for the function x defined in Theorem 4.1. We will have use for the fact that 
X has an elementary antiderivative x^ an d a double antiderivative x^ given 

by 

a 1 smh(7rx/h ) 

X ( x ) = o„u ^ 

(4.4) 


87 rh cosh( 7 Tx/h) + cos( 7 tQ) ' 


1 


X^{x) = -3- log(cOSh(7rx//l) + COs(7T0)), 

87T z 

respectively. While there in general seems to be no nice closed form of the 
leading order surface profile 

m = -(9 - a 2 ^) _1 x 


obtained in Theorem 4.1, we can still give some of its properties. An immedi¬ 
ate one is that 772 is smooth. In Proposition 4.4 we give a series expansion 
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for y 2 in powers of e 1T \ x \/ h . Furthermore, perhaps more surprisingly, we can 
find an explicit expression for 772 in terms of elementary functions whenever 

(4.5) m := 

TTQL 

is a natural number. If m £ N, then e^^ x l a = e ±m7rx / /l are integral powers 
of e :iz7TX / h , which would appear on the right side of Equation (4.3) if we had 
written out cosh(7r x/h) and sinh(7r.T/h). Since x 1 —> spans the kernel 

of g — a 2 d 2 , this explains integral values of rri being special. 

Before we state Proposition 4.4 and Theorem 4.6, we need a lemma to 
simplify some expressions. 

Lemma 4.3. For m £ (0, 00 ) \ N and 9 £ (0,1), we have 



(4.6) 


1 OO 

— + 2m V (—1) 

m 


k cos(kir9) cos(m,ir9) 

= 7T- 


fc=l 


m 


2 _ k 2 


sin(m7r) 


which is equal to 
(4.7) 


[°° v m- 1 cos(ir9)y + 1 
/ 0 y 2 + 2cos(7T0)y + 1 


whenever m £ (0,1). Fwdhermore, for m £ N, 


fc=i 
k^m 


. , COs(m7T0) „ . , 

-(—1) -- + tt 6 sm(m7T0) 

2m 


Proof. Both sides of Equation (4.6) define meromorphic functions in m on C 
with simple poles in the points Z x {0}. Moreover, they are both equal to the 
integral in Equation (4.7) when m £ (0,1), which can be seen by calculating 
the integral with both the residue theorem (around a keyhole contour) and 
a Laurent series expansion of the integrand. Since the interval consists of 
non-isolated points, we have equality on all of C. Finally, Equation (4.8) 
follows from Equation (4.6) by taking limits. □ 

Proposition 4.4 (Expansion for 7 / 2 ). If the number m in Equation (4.5) 
satisfies m £ (0, 00 ) \ N, then the leading order term . of the surface profile 
from Theorem f.l is given by 


m(x) = 


1 


87r 2 a 2 


log(l + 2cos(7r0)e-" |a;|//l + e~ 2n W /h ) 


_ T Cos(m^) e _ v5 |,| /a + 2m2 y t COS(farf) 

1 m / 7 r 1 ' h(m^ — \ 


sm(m7r 


k =1 


k(m 2 — k 2 ) 
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while if m £ N, then 

1 


mi x ) = 


8ir 2 a 2 
+ (-!)' 


log(l + 2 cos ( t T6)e-*W /h + e~ 2 *W h ) 
3 cos(imrd) 


7T\X\ 


2 771 


+ 7T0 sin(m 7 T 0 ) + cos(m7T0)^ ) e~^\Ah 




fc=l 


These series converge uniformly, and, excluding the origin, so do the series 
for the first, derivative. Moreover, when m £ N, the function y 2 is given 
explicitly in terms of elementary functions by 


mix) = 


1 


87 r 2 a 2 


m—1 


1 _ k cos((m — k)ir6) , 

+ 2 2 ^(-l) m —-—- cosh((m - k)irx/h) 

k= 1 


m 


+ r(e 7rx ^ h ) + r(e ~ nx/h ) 


is defined by 


where r: ( 0 , 00 ) 

r(x) := ^( — l) m cos(mir9)x~ m log(l + 2 cos( 7 T 0 )x + ,x 2 ) 

+ (—l) m sin(m7r0)x~ m (arctan(cot(7r#) + csc(7T0)x) — 7r( 1/2 — 6)). 
Proof. It follows from 

'2 a 


V 7T ( 


,2 1 £2 ! a > 0) 


■ a 2 + £ 2 

and the definition of 772 , that we may write 772 as the convolution 


mix) = -v-*-=(e v^M / 0 * x) 0 ) 

2 cxyfg 

1 

iJ{x,x) + J(~x,x)), 


(4.9) 
where 

Equivalently 

(4.10) 

(4.11) 


2 ay/g 

J(x, X )--=e-^ x,a f X e^v/ a X {y)dy. 


mix) = i x i *) + J i~ x ’ x 

= ~2X^( x ) - 7^1 + ^(-^x" 


a*'~ ' ' 2a 3 

through integration by parts, where and are the antiderivatives defined 
in Equation (4.4). 

We first use Equation (4.10) to obtain an explicit expression for 772 when 
m £ N. By using the substitution x i->- e nx ^ h , we find that 

Ax,x") = Afiie 7rx/h ), flix) := /V'-iiTTirn 

8ir z J o z z + 2cos(7rtl)2; + 1 


dz. 
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The fraction in the integrand in the definition of f± has partial fraction 
decomposition 


z 2 -l 


AttQ 




z 2 + 2cos(7 t6)z + 1 


= 1 - 


z + e 


inO 


z + e 


—inO ’ 


and since 


1 


z + a 

this means that 


(— a) m 
z + a 


m —2 

£ 

k =0 


(-a) 


m—k—1 


z K , aGC^ / —a, 


/i(x) = — \-(—l) m e imnU x~ rn \og(x+e mU )+(—l) m e~ im7rU x~ m \og{x+e~ mt> ) 
m 


g—imnO ^—m 


fc=l 


where log(-) denotes the principal branch of the logarithm. The result now 
follows by using the identity 


log (a; + e™ 6 ) = ^ log(l + 2 cos(7r#)x + x 2 ) 

— i(arctan(cot(7T0) + csc(7T0)x) — n/2), 


valid for all x£l. 

For the series representation of r/ 2 , we use Equation (4.11), because this 
leads to a series that converges more rapidly. We will assume that m £ 
(0, oo) \ N; the case for rn £ N is similar, except that one needs to use 
Equation (4.8) instead of Equation (4.6). We use the same substitution as 
before to arrive at 

where fa'. (0, oo) —> M is defined by 

f 2 (x) := mx~ m f z m ~ l log((,z -1 + z)/2 + cos(7T0)) dz. 

Jo 

One may check that one has 

log((2 -1 + z)/2 + cos(t t6)) = - log(2) - log(z) - 2 ^(-l) fc COS ^ 7r6> ^ fc 

k =1 

for z £ (0,1) and 

logtC*- 1 + z)/2 + cos(vr0)) = - log(2) + log(z) - 2V(-1 z ~ k 

k 

k =1 


for z £ (1, oo). 

It then follows by termwise integration that 


f2{x) 


1 

m 


log(2) - log(x) - 2m^(-l) fc 
k =1 


COS(/C7T0) k 

k(m + k ) 
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on (0,1] (the endpoint is Abel’s theorem [34, Theorem 17.14]), and that 
/ 2 (x) = h(l)x~ m + mx~ m J* z m ~ l log((z + z - 1 )/2 + cos(t t 6)) dz 

= - — - log(2) + log(x) - 2m^(-l) fc cos ^’ 7r ^ 

m z J 


k =1 
oo 


m 


+ 4m »1) 


k =1 


k(m — k) 

k cos(knd) 
m 2 — k 2 


x 


for x € [l,oo). Employing Equation (4.11), we find that 772 is given by 

1 


m(x) = 


8ir 2 a 2 


log (cosh (irx/h) + cos(7r#)) — 7r|x| /h + log(2) 


m 


+ 2m^(-i) 


_ n fc cos(fc7r0) \ -yglxl/a 


fc=i 

+ 


m. 


2 — A; 2 


2m 2 V(-l) fc , roS 2 (fc7r ^ 

' himA — 




fc(m 2 — k 2 ) 


for all x 6 1, by using that 772 is even and observing that for x > 0 we have 
gTrrr/h £ [l,oo) and e~' KX l h £ (0,1]. If we now apply Equation (4.6) from 
Lemma 4.3 in order to get a closed-form expression for the coefficient in front 
of e^v^M/ a ; we arrive at the desired expansion. □ 


Remark 4.5. The only obstacle to convergence of the series given in Propo¬ 
sition 4.4 is the origin; thanks to the exponential factor e ~ kir \ x \/ h ^ the con¬ 
vergence is rapid away from the origin. It should also be noted that, while 
Equation (4.3) seems to suggest that 7/2 should be expandable in a series 
of powers of sech(7rx//i) by equating coefficients in the differential equation 
defining it, this seems to lead to a series that does not converge. We have 
kept the series expansion for 772 also when m £ N, because the expression 
in terms of elementary functions is unwieldy, and prone to numerical errors 
even for small values of m. 


The expressions found in Proposition 4.4 have well defined pointwise limits 
as 9 f 1 (for r/0) and 0 4-0- In particular, when m. = 1 these are given by 

UmT72(x) = —^2 [1 - e^logU + ^ x/h ) - e~™/ h log(l + e™/ h )\ 

lini77 2 (x) = —^[1 + e KX / h log 11 - e~™/ h \ + e~™' h log 11 - e™/% 

6»tl 87T-Q;- 

which can can be seen as graphs drawn with thicker lines in Figure 1, together 
with 772 for various values of the parameter 0. 

We see from Figure 1 that one gets a depression at the origin, which 
becomes more pronounced the closer the point vortex is situated to the 
surface. The profile when the point vortex is close to the surface is very 
similar to the profile for the infinite depth case, found in [41]. However, 
a feature which is not seen on infinite depth is that there is a significant 
difference between the case 6 < 1/2 and the case 6 > 1/2 (in addition to the 
changing sign of ci). For 8 < 1/2 there is a single trough at the origin, and 
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FIGURE 1. The leading order term 77 2 in 77 (e), with h = 
l,a 2 = 1 /( 87 r 2 ),m = 1. The values of 9 shown are 9 = 
0.1,0.2,..., 0.9, together with the thicker lower and upper 
limits 9 4 . 0 and 9 f 1. 


772 is everywhere strictly negative. When 9 > 1/2 one in addition gets crests 
on either side of the origin. As we can see from Figure 1, the positions of 
these crests depend on the position of the point vortex. 

Some of what we have just discussed is not limited to the specific choice of 
constants that are used in Figure 1, and for which Proposition 4.4 yields an 
explicit expression for 772 . We will see that m = 1 plays a special role in the 
asymptotic behavior of 772 , however. More precisely, we have the following 
theorem: 


Theorem 4.6 (Properties of 772). The leading order surface term 772 always 
satisfies 772(0) < 0 and 772(0) > 0, meaning that the origin is a depression. 
When 9 < 1 / 2 , the function 772 is everywhere negative, and strictly increasing 
on [0, 00 ). For 9 > 1 / 2 , we have two cases, depending on the number m 
defined in Equation (4.5): 

(i) If 777 > 1 /( 2 9 ), then 172(3:) is positive for sufficiently large |x|. In 
particular, 772 has crests on either side of the origin. 

(ii) If m < 1 /( 2 9 ), then 772(37) is negative for sufficiently large |x|. 
FuHhermore, 772 has the following asymptotic properties for any 9 £ (0, 1 ): 

(i) For 7 ? 7 . > 1 


(4.12) 


lim i)2{x)e KX ^ h = 

x—foe 


2 COs( 7 T$) 
m 2 — 1 87 r 2 a 2 


(ii) If m = 1, then 

e"/ ft COS( 7 T 0 ) 

Inn 772 ( 3 ;)-— = — , 9 . 

»oo -Kx/h 87r-or 


lim m {x)e^ x / a = 


7r cos(mn9) 

sin(?777r) 87T 2 o 2 


(Hi) For 777. < 1 
(4.13) 
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Proof. We first prove that 7 / 2 ( 0 ) < 0 and 7 / 2 ( 0 ) > 0, which holds for all values 
of 77 ;. and 9. By inserting x = 0 in Equation (4.9), and using the evenness of 
X, we find 


772 ( 0 ) 


1 f°° 

—-/ e-^y/ a x (y)dy 
Uy/g Jo 

- , / e-^y/ a x\y) dy < 0 , 

a Jo '-v-' 

> 0 on (0, 00 ) 


where the second equality follows from integration by parts, and the function 
X^ was dehned in Equation (4.4). Since 7/2 = —(g — a 2 d^.)^ 1 x, we also have 

772 ( 0 ) = ~2 ( 5772 ( 0 ) + x( 0 )) 

or 

= 2 a 3 j e _ ^ |y|/a (x(°) ~x(y))dy 
> 0 , 


as x achieves a global maximum at the origin. 

Suppose now that 6 < 1/2. Like in Proposition 4.4, we use the fact that 
7/2 may be written as the convolution 

(4.14) ?/2 = — 2ax^ e_v/ ^ 1/0 * X ' ) ’ 

which shows that 7/2 is strictly negative, since x is strictly positive when 
9 < 1/2. Moreover, some manipulations of the above formula shows that we 
may write the derivative of 7/2 as 


7?2 (" 7-0 = - 


ay/g L 


sinh J e '^ )y ^°"x(y) dy 

+e -pgx/a J ginh x\y) dy 


where we have used the fact that % is even. One may check that x! is strictly 
negative for x > 0 when 9 < 1/2. This shows that i/ 2 is strictly positive for 
x > 0, and so 7/2 is strictly increasing on [0, 00) by the mean value theorem. 

Before we consider the case 9 > 1/2, we prove the asymptotic properties 
for 7/2 listed in Equations (4.12) to (4.13). These follow by multiplying each 
side in Equation (4.9) with the appropriate factor and taking limits. For 
instance, suppose that 777 > 1, meaning that y/g/a > 7 r/h. For the integral in 


g nx/h 



e V9v/ a x{y) dy 


fl oo e'/y yla x{y) dy 

e {y/d/oi—n/h)x 


there are two possibilities: If 9 = 1/2, then it is possible that the integrand 
is integrable on the entire real line, meaning that the limit as x —)• 00 is zero; 
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otherwise, the integral tends to Too, and so 

v I- 00 e ^ y/a x(y)dy e^ x/a x(x) 

lim -— lim - 

x —>-oo (X\/9/ Q —^/E)x x ^ \J~y/Oi — 11 / hf) ^ a 7r/h)x 

1 COs(7T0) 

^/g/a-ir/h Ah 2 

by L’Hopital’s rule. The other limits can be treated in a similar way, with 
one exception: 

The procedure will show that when m < 1, we have 

1 roo 

lim rj 2 (x)e^ x ^ a =-/ e^ y ^ a x(y) dy 

x ^°° 2 aty/g J _ 00 

1 f°° m cos(irO)y 2 + 2y + cos(7T0) 

87 T 2 a 2 m J 0 ^ ( y 2 + 2 cos (7 xd)y + l ) 2 ^ 

= _ 1 rV-l CQS(7T% + 1 

87r 2 a 2 J o y 2 + 2 cos(7t0)t/ + 1 

where the second and third equality follows from the substitution y 1 —> e wy ^ h 
and an integration by parts, respectively. The result now follows since the 
integral on the final line is equal to the right-hand side of Equation (4.6) by 
Lemma 4.3. 

Finally, we consider the case of 9 > 1/2, which is harder to describe 
completely, as the integrand in the convolution in Equation (4.14) changes 
sign. Observe that the claims on the sign of 7 / 2 ( 2 :) for sufficiently large x 
follows for m j- 1/(20) from the limits in Equations (4.12) to (4.13). An 
additional argument is needed for the edge case ru. = 1/(20), because the 
limit in Equation (4.13) vanishes. It turns out that Equation (4.12) also holds 
in the special case m = 1/(20), which can be shown with the same method 
we used to show the other limits. Hence 7/2 is negative for sufficiently large x 
when m = 1/(20), which exhausts the values of m. □ 

Remark 4.7. It is likely that 7/2 has similar properties to those for the case 
0 < 1/2 when 0 > 1/2 and m < 1/(20), but we have not been able to prove 
this. 

We are now in a position where we can give the sign of C 3 in the expansion 
in Theorem 4.1 for 0 < 1/2. 

Proposition 4.8 (Sign of C 3 ). The constant C 3 in Equation (4.1) is negative 
when 0 < 1/2. In particular, if 0 = 1/2 and e is sufficiently small, the waves 
obtained in Theorem f.l are left-moving when e > 0 and right-moving when 
£ < 0 . 

Proof. Recall the definition of ((3 in Equation (4.1). From Theorem 4.6 we 
know that 7/2 is negative, and strictly increasing on [0, 00 ). Furthermore, the 
factor ci + 4>j/(-, 0) is positive and strictly decreasing on the same interval. It 
follows that also £3 is positive and strictly decreasing on [ 0 , 00 ). 

The harmonic function H{ 0)^3 on M x (— h, 0) assumes the value 0 at 
the bottom of the domain and £3 > 0 at the top of the domain. By the 
maximum principle, it is positive on the entire domain. Thus we may use 
the Hopf boundary point lemma (see [23, Lemma 3.4]) in order to conclude 
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that [/f(0)C3]j/(0, — h ) > 0. The result will therefore follow if we can show 
that [H(0)^] y is increasing along the y-axis. We will do this by looking 
at [H(0)(^3] x on (0,oo) x (— h, 0). Because of its values on the boundary, 
it is negative in the interior. Another application of the Hopf boundary 
point lemma implies that [i?( 0 )^ 3 ] xx is negative on the y-axis (except at the 
point (0 ,—h), where it vanishes). Since [H(0)C^] yy = — [if(0)£3] xx by the 
harmonicity of H( 0)^3, this concludes the proof. □ 

We finish our exposition on a single point vortex with a short discussion on 
the streamlines of waves obtained in Theorem 4.1. Observe that if (. x(t ), y(t)) 
denotes the position of a fluid particle at time t, then 

(4.15) (x(t),y(t)) = w(x(t),y(t),t), 

before the new variables in Section 2.2. After introducing the steady variables, 
Equation (4.15) becomes 

(4.16) (x(t),y(t)) = w(x(t),y(t )) - (c, 0), 

meaning that if we only keep the first order terms for w and c from Theo¬ 
rem 4.1, we obtain (keeping the same notation for the paths) 

(4.17) (x(t),y(t)) = eV 1 - ($ + ay) (x(t),y(t)). 




(a) e = 1/3 (b) e = 2/3 

Figure 2 . Streamlines in the frame of reference traveling with 
the wave, for h = n and e > 0. The wave corresponding to 
8 = 1/3 propagates to the right, while the wave corresponding 
to 9 = 2/3 propagates to the left. The arrows illustrating 
the vector held on the right hand side of Equation (4.17) 
have been scaled here for visibility, and only their direction is 
quantitatively accurate. 


We have used this to obtain Figure 2, which shows streamlines in the 
steady frame moving with the wave. The portraits corresponding to 9 and 
1 — 6 can be obtained from each other by a 180° rotation. When 9 = 1/2, all 
the streamlines are closed (not shown), so we will focus on the case 9 ^ 1/2. 
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The lines y = —h and y = 0 are nullclines for the system in Equation (4.17), 
and the points (x. y) with 

x = ±h/7r arcosh(|2sin(7T0) tan(7rd) + cos(7rd)|), 

(4.18) _ i-h 9 < 1/2 

V ~ [0 9 > 1/2 

are equilibrium points, corresponding to stagnation points. One may check 
that 

h/ir arcosh(2 sin(7T0) tan(7rd) + cos(7T0)) = \/3h9 + 0(9 5 ) 

as 9 4- 0, meaning that the distance between the equilibria is very close to 
linear in 9 for small 9 (a corresponding statement holds for 1 — 9 small). They 
go off to infinity as 9 —> 1/2 from either side. The heteroclinic orbit (which 
can be expressed explicitly in terms of arcosh) connecting the two equilibrium 
points described in Equation (4.18) encloses a critical layer containing closed 
streamlines. Outside this region the particles always move in the same 
direction with respect to the steady frame. This direction is either to the left 
or right depending on the sign of cot(7rd) and e. 

We also mention that on infinite depth, the streamlines always look like 
those in Figure 2b. If the point vortex is situated at (0, — d), the equilibrium 
points at the surface will be at (±\/3d, 0), and the points on the heteroclinic 
orbit between these satisfy 

x 2 + (y + d) 2 = 2dy(l + coth (y/ (2d))), 

which is close to half an ellipse centered at (0, —d) with semiaxes \/3 d and 
~ 2.0873d. The equilibrium points in Equation (4.18) converge to those on 
infinite depth as h —> oo if d is held fixed. 

Because only the first order terms in e have been kept in Equation (4.17), 
we do not make any claim about the accuracy of the phase portraits in 
Figure 2 for the full system in Equation (4.16). That would require further 
and more thorough analysis, in particular for the case 8 = 1/2. Still, the 
phase portraits can give some indication as to how these waves look beneath 
the surface. One feature will remain the same for Equation (4.16): Because 
of the singularity of <f> at (0, —(1 — 9)h), the streamlines will always remain 
closed sufficiently close to the point vortex. 

5. Several point vortices 

We aim to extend the existence result for traveling waves with a single 
point vortex in Theorem 4.1 to a finite number of point vortices on the y-axis. 
As opposed to the single vortex case, where we could choose 9 freely, there 
will be limitations on the positions that the point vortices can occupy. We 
will return to this. Suppose that 

1 > 6\ > 6 2 > ■ ■ ■ > 0 n > 0 , 

and that we wish to establish the existence of a traveling wave with point 
vortex at the points 


(0, _(l - dr)/,),..., (0,-(!-£„>) 
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the situation being otherwise similar to that of a single point vortex. The 
admissible surface profiles are those in A q x , as the uppermost point vortex is 
the most restrictive. 

For r\ £ and 7 = ( 7 1 ,..., 7 ”) £ R n we may define 

n 

(5.1) <|F : ^V<|£. 

j = 1 

where 

&( X y) := JL log ( '^M ~ , 7 = i ... n 

’ 47t \cosh(7rx//i) + cos(ir(y/h + ftj)) / ’ ’ 

in f2(r/). We will seek solutions of the form 

w = V ± [H{r))( + &], 


cf. Equation (3.6) for a single point vortex. 

The main difference from the single point vortex case is of course the 
vorticity equation, Equation (2.9), which needs to be imposed for each of the 
point vortices. For the ?'th point vortex, the vorticity equation reduces to 


M) 


= V ± [H( V )(](0,-(l-9 i )h) 

1 / n 

+ 4ft ( 7 * cot ( 7r6l i) + 

k i =1 


cot 





which, if we assume that i] and are even (see the discussion before Equa¬ 
tion (3.7)), can be written more succinctly as 


(5.2) cl = -([H(r,)() y (0, -(1 - 9i)h))? = 1 + 0 7 . 


Here, we have defined 1 := (!,...,!) £ R n and the matrix 0 £ R nxn by 


(5.3) 


&i,j — \ 


4 h 
1 


COt(7T0,;) 


- , + , &i + 0j\ ( &i — Oj 

4ft- V V 2 ) I 2 


i = j 

i + j 


for 1 < i,j < n. 

As opposed to for one vortex, it is now more natural to use the wave velocity 
c as the bifurcation parameter. We will therefore write e instead of c in order 
to have notation that is more consistent with the one vortex case. The idea 
is to use the vortex strengths 7 in order to balance Equation (5.2), which is 
possible when 0 is invertible. It should be emphasized that this is almost 
always the case (Theorem 5.6), but that there always are configurations of n 
point vortices that yield singular 0 (Proposition 5.7). We have already seen 
such a configuration, albeit a trivial one: For the case n = 1 one has 0 = 0 
when 9 = 1 / 2 . 

We make the necessary redefinitions 


X s := H s even ( R) x H s even ( R) x R ? \ 
W := H!~l( R) x H s even ( R) x R", 
U$ i :={(r,,(,l)£X s :r 1 eAe 1 } 1 
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and proceed to define, for s > 3/2, the map F\ : Ug xl-) i?| ve ^(M) by 


= e 


v'C + G(r,K 


+ <f > 7 

W ) 2 y \ 

+ (C' + (i,,/).v^ + ( 0(; ,K + (-/.i)-v^ + _ )f 


2W> 


the map F 2 : xRa i7| ven (M) by 

^ 2 ( 77 , C, 7>e) := £h + C + 4> 7 , 


and finally the map T 3 : x M —x M n by 

F 3 {v, C> 7, e) : = ©7 - el - ([#fa)C]i/(0, -(1 - 0i)/O)"=i- 
In all of these dehnitions, the function <J > 7 and its derivatives are evaluated 
at (x,r/(x)), which is suppressed for readability. 

We now define F := (T\, F 2 , T 3 ): Ug —X Y s , and seek solutions of the 
equation 


(5.4) 


^(^,C,7,e) = °- 


which has the origin as a trivial solution. We are led to the following analog 
of Theorem 4.1 for several point vortices, establishing the existence of a 
family of small, localized solutions, assuming that 0 is nonsingular. The 
resulting waves have one critical layer for each point vortex, assuming that 
no component of 7 vanishes. 


Theorem 5.1 (Traveling waves with several point vortices). Let s > 3/2, 
and let 1 > 9\ > 62 > ■ ■ ■ > 9 n > 0. Suppose that the matrix 0 defined in 
Equation (5.3) is invertible. Then there exists an open interval I 3 0 and a 
C°°-curve 

I -x (H s even (R) n A*) x H s even (R) xrxl 

e 1 x (r)(e), C00,700,<0 

of solutions with velocity c = e to the Zakharov-Craig-Sulem. form.ulation, 
Equation (5.4), for point vortices of strengths 7 1 (e),... ,7 n (e) situated at 

(o,-(i-^)h)„..,(o,-(i -e n )h). 

The solutions fulfil 

77 (e) = 772 e 2 + 0 (e 4 ), 

(5.5) C( e ) = ( 3 £ 3 + 0(e 4 ), 

7(e) = 71 e + 7 3 e 3 + 0(e 4 ), 

in their respective spaces as e —X 0 , where 71 := 0 _ 1 1 , the function 772 6 
i7| wen (M) is defined by 

m ■= -(ff-« 2 ^) _ 1 X, X := 4> 71 (-,0) + ^ 71 (-,0) 2 , 

and where 


C 3 = - ? 72(1 + 4> 71 (.,0)), 

73 = 0 _1 ([-f^(O)C 3 ]y(O, —(1 — di)h))2 =1 . Y 
with <h 71 as in Equation (5.1) and H as in Definition 2.2. 
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Moreover, there is a neighborhood of the origin in Uq X K such that this 
curve describes all solutions to F(rj, £, 7 ,e) = 0 in that neighborhood. 


Proof. As for a single point vortex, we wish to apply the implicit function 
theorem at the origin. We find the derivative 


D X F{ 0,0, 0,0) 


g- a 


2Q2 


0 

0 


0 


J W S ve„(R) 

-([tf(0)-],(0,-(l-^))?=i 


0 

0 

0 


where ([i4(0)-] y (0, —(1 — 0 ,:)/i ))" =1 means the operator i4| ven (M) —> M n defined 

by 

C ([H(0)C) y (0, -(1 - 9i)h))^ =1 . 

Recalling that g — a 2 and 0 are invertible by the discussion after Equa¬ 
tion (4.2) and by assumption, respectively, DxF{ 0, 0, 0,0) is an isomorphism. 

Hence we can use the implicit function theorem to deduce the existence of 
an open interval / around zero, an open set V C Uf containing the origin, 
and a map / G C°°(I, V) such that for ( rj , f, 7 , e) G V X /, we have 


F(r], C, 7 , e) = 0 (77, C, 7) = /(e)- 


The terms in the expansion in Equation (5.5) can be obtained as in the proof 
of Theorem 4.1. □ 


Remark 5.2. It is worth mentioning that on infinite depth, the matrix 0 is 
always invertible. A corresponding existence theorem for infinite depth would 
thus hold for any configuration. 

Rem.ark 5.3. An extension of the existence result in Theorem 5.1 to point 
vortices that are not all on the same vertical line would require a different 
argument than the one we have used. The main issue is that assuming 77 and 
( to be even is then no longer sufficient to satisfy the vertical component of 
the vorticity equation, like we did to obtain Equation (5.2). 

One may note that the sign reversal of the wave velocity about the midpoint 
6 = 1/2 that we saw with the single point vortex, can be seen also for several 
point vortices, albeit in a different manner. If the matrix 0 corresponds to 
1 > 6\ > • • • > 6 n > 0, and we reflect the vortices across the line y = —h/2 
by considering di := 1 — 6i, 1 < i < n instead (without reordering them), 
then the new matrix is —0. This causes a swap of sign on the leading order 
vortex strengths, 71 = © - 1 1 . 

We have pointed out that the matrix 0 is not invertible for all configurations 
of point vortices, and gave the trivial example of 6 = 1/2 for a single point 
vortex. This example, together with Theorem 4.1, also shows that invertibility 
of 0 is not a necessary condition for the existence of a traveling wave with 
point vortices in those points. See also Remark 5.5. 

The only case for multiple point vortices on the y-axis where we can 
feasibly describe the admissible positions directly is for n = 2. In fact, we 
give a complete description of when 0 is invertible in Proposition 5.4; see 
also Figure 3, which presents this result graphically. One may observe that 
the midpoint between the bottom and surface plays a role also here. 

Proposition 5.4 (0 for n = 2). For two point vortices, we have the following: 
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Figure 3. The determinant of 0 for the case n = 2 as a 
function of (9 1 ,^ 2 )- The determinant vanishes along the solid 
black curve, which is given explicitly as a parametrization in 
Proposition 5.4. (In the figure, the level curve for det(@) = 0 
is computed numerically.) 


(i) If 6 \ < 1/2, then 0 is invertible for all 62 £ (0, 9\). 

(ii) If 9 1 > 1/2, then © is invertible for all 62 £ (0, 9\) except for exactly 
one value, 0 < ^ 2 (^ 1 ) < 1/2. The graph of 9 2 : (1/2,1) -A (0,1/2) is 
described by the curve 

(tt/4, 3tt/4) -a (1/2,1) x (0,1/2) 

t ^ (f + /(f),f-/(f))/7T 

where f: (71 /4,37T/4) -A M is defined by 

f(x) := arccot (^ ^cot(x) 2 + y/4 — 3 cot(x) 4 ^ \ . 


Proof. It is useful to write the determinant of 0 as 


det(0) 


1 

16 h 2 


COt(7T0i) COt(7T02) + 


4 sin(7T0i) sin(7T02) 

(cOs(7t 6*2) — COs(7T0i)) 2 


One immediately observes that the second term inside the parentheses is 
always strictly positive. If 9\ < 1/2, then one has in addition that the first 
term is nonnegative for any 62 £ (0,#i) Q (0,1/2). This proves the first part 
of the proposition. 

For the second part, let us first prove that there is exactly one value of $2 
for each 9\ £ (1/2,1) that makes 0 singular, and that this value lies in the 
interval (0,1/2). For fixed 9\ £ (1/2,1) the determinant is strictly increasing 
in 62 , and tends to —00 as $2 i 0, and to 00 as $2 t $i- Hence it vanishes at 
exactly one value of 62 , say 6 * 2 (0i )■ Because the determinant is positive when 
$2 = 1 / 2 , this value must necessarily he in the interval ( 0 , 1 / 2 ). 
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We now move to the parametrization of the graph of the map 62 '■ (1/2,1) —> 
(0,1/2). One may note from Figure 3 that there is symmetry across the 
diagonal line 

{(01,02) £ (0, l) 2 : 0i > 02, 01 + 02 = 1}, 

which suggests making a change of variables. By letting 


(5.6) 


4>1 : = 7T 


01 + 02 


4 >2 : = vr 


01 - 02 


we can write the determinant in the form 


det(0) 


1 

IQh? 


COt(</>i) 2 COt(</>2) 2 — 1 
COt(</> 2 ) 2 - cot(0i) 2 


+ cot( 02) 2 


COt(0l) 2 , 


which leads us to solve the quadratic equation 

x 2 — ax + a 2 — 1 = 0, a := cot 2 (0i), x := cot 2 ( 02 ) 


for x, given a. Doing this yields the parametrization, by using <f>\ as the 
parameter (some care has to be taken to ensure that one picks the right 
branches of the functions involved) and going back to the original variables 
by inverting Equation (5.6). □ 


Remark 5.5. By employing the parametrization of the graph of 02 provided 
by Proposition 5.4, one can show that the each column of 0 is linearly 
independent from 1 when det(0) = 0. This implies that an argument similar 
to that of Theorem 5.1 can be performed, by using the vortex strength 7 1 as 
the bifurcation parameter, instead of c. Thus it is possible to show existence 
for any configuration when n = 2. An extension of this argument to n > 2 is 
harder, because it requires the rank of 0 to be n — 1 . 


While the set of configurations that make det(0) vanish is hard to describe 
in general when n > 2, some observations can be made. Of course, if n > 2, 
and as long as the derivative of det( 0 ) with respect to the variable ( 0 i,..., 9 n ) 
does not vanish at a point where det( 0 ) = 0 , the zero set of det( 0 ) is locally 
a smooth manifold of dimension n — 1 around that point by the implicit 
function theorem. When n = 2, the zero set is actually the graph of a smooth 
function in 6 \ by Proposition 5.4, and numerical evidence suggests that the 
zero set is the graph of a smooth function in (0i, 62 ) when n = 3. Actually 
checking that the derivative does not vanish is hard, but we have the following 
theorem: 


Theorem 5.6. The subset of configurations of point vortices in 
{(01,..., 0 n ) G (0, l) n : 1 > 0i > 0 2 > • • • > 9 n > 0} 
such that 0 is not invertible has measure zero. 


Proof. Each entry in 0 is analytic in each 0j for 0i,..., 0j_i, 0,:+i,..., 9 n 
fixed. It follows that det(0) also has this property, when viewed as a function 

U := {(0i,... ,0 n ) G ( 0 , l) n : 1 > 0i > 0 2 > • • • > 0„ > 0 } -> M. 

We first verify that det(0) does not vanish identically on U. To that 
end, fix 1/2 > 0i > 02 > • • • > 0 n > 0 and consider 0i = e9\,... 9 n = e9 n 
for 1 > £ > 0. The purpose of the upper bound of 1/2 is to make sure 
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that tan(vr0j) is well defined for all 1 < i < n. Observe now that if we let 
T := diag(tan( 7 T 0 fc))£ =1 , then 


4h[T0 } itj = 


where 


limtan(e7T0j) cot sn 


f 1 

|tan( 7 T 0 j) (cot (4— cot ^ 7 r 

i ( e i + §; ' ' ■ 


di—dj 


l = J 

i 7 k j 


— cot STT 


6i - 6 j 






for i j. It follows that diag(tan( 7 T 0 fc ))£ =1 0 has a limit in BlW 1 ) as e 4- 0, 
and that this limit is 


(5.7) 


limdiag(tan( 7 r 0 fc ))fc =1 © = - B), 

£4-0 4 ri 


where we have defined B £ M nxn by 


(5.8) 


Bij 


0 

Midjift - e])- 1 


i = j 

i + J 


In particular, B is skew-symmetric, which implies that — B is in¬ 
vertible. Since the set of invertible operators is open, so is the matrix 
diag(tan(7T0fc)))! =1 0 for sufficiently small e. which in turn means that 0 is 
invertible for such e. 

Finally, the set U is connected. Hence, since we know that det(0) is 
analytic in each variable and does not vanish identically, we infer 5 that the 
subset of U on which det(0) vanishes has measure zero. □ 


In general we cannot do better than Theorem 5.6, in the sense that for 
any n > 1 there will always be a configuration of n point vortices that makes 
det( 0 ) vanish. 

Proposition 5.7. There are always configurations of point vortices in 

{( 01 , • • •, 0„) £ ( 0 , l) n : 1 > 0 i > 6 2 > ■ ■ ■ > e„ > 0 } 

where 0 is singular. 

Proof. The matrix appearing on the right-hand side of Equation (5.7) in 
the proof of Theorem 5.6 has a positive determinant. Indeed, the matrix 
B defined in Equation (5.8) is skew-symmetric, so its spectrum is purely 
imaginary. Moreover, since the matrix is real, the eigenvalues are either zero 
or appear in complex conjugate pairs. 

Say that the first m eigenvalues Ai,..., A m of B are zero and that 

Ato+2j — 1 — ^m+2j — iTji j — ^l)/2, 

where the fij are real. Then it follows that 

det (-^(/r» - B)^j = det(/ Rn - B) 

= ( 4 / l )n ( 1 + + M2) •••(!+ M(n- m )/2)1 

5 This follows by induction on the dimension, by using the well known result in one 
dimension. 
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because the determinant of a matrix is equal to the product of its eigenvalues 
(taking algebraic multiplicity into account). By Equation (5.7) we then have 

n 

(5.9) det(0) tan(7T0fc) > 0 

k =l 

for small e > 0 (as in the proof of Theorem 5.6) by continuity of the 
determinant. Since all the tangents are also positive, this implies that 
det(0) > 0 for small e > 0. 

It remains to exhibit a configuration where det(0) < 0. To that end, fix 
^ > 8 \ > §2 > O 3 > ■ ■ ■ > 0 n > 0 and consider 0 \ = 1 — e 0 \. O 2 = £ 62 , • •. 0 n = 
e0 n for 1 > e > 0. Proceeding as in the proof of Theorem 5.6 we find 

(5.10) limdiag(tan(7T0fc))£ =1 0 = - B), 

£ 4,0 4 11 


where we have defined B E R nxn by 


Bi,j ■= 



i = j or i = 1 or j = 1 , 
0|) _1 otherwise. 


This matrix is still skew-symmetric like B, and so the right-hand side of 
Equation (5.10) has a positive determinant, as before. Hence Equation (5.9) 
holds for small e. However, now tan(7T0i) is negative and the rest of the 
tangents are positive, meaning that det(0) must be negative. □ 


6. Explicit expressions for infinite depth 


In this section we give some explicit expressions for periodic waves with 
a point vortex on infinite depth, constructed in [41]. We will adopt the 
notation and conventions used there. The fluid domain for the trivial surface 
is R X (— 00 ,1) and the waves have period 2nL. The stream function for the 
rotational part is denoted by G. 


Proposition 6.1 (Stream function). The stream, function for the rotational 
part is given by 


G{x,y) 


1 f cos(x/L) — cosh(y/L) \ 
47 t g ycos(x/L) — cosh((y — 2)/L)J 


Proof. We wish to find the stream function G:lx (— 00 ,1) —X M corre¬ 
sponding to equally spaced point vortices of unit strength at the points 
2'kLT, X {0} C R 2 , and which is such that this stream function vanishes at 
the surface, R X {1}. By symmetry, it must be the case that G x vanishes on 
7 tL(1 + 2Z) X (— 00 ,1). This leads us to the boundary value problem 

AG = 5, G|^ =0 = 0, G x \ x=±7rL = 0, 

on (—7 tL, nL) X (— 00 ,1). This equation can be dealt with using Theorem A.2 
in Appendix A. 

In order to apply Theorem A.2 we require a conformal map satisfying the 
requirements in the theorem statement. One may check (see [45, Sections 7.1 
and 7.2]) that 

._ tanh(l/(2L)) — tanh((l + iz)j (2L)) 
tanh(l/(2L)) + tanh((l + iz)/(2L)) 
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defines a bijective conformal map from the half strip (—ttL,ttL) X (— 00 ,1) 
onto the slit unit disk D \ ((0, exp( — 1/L)) X {0}), and which is such that 

(i) The origin is fixed. 

(ii) The surface is mapped to the unit circle. 

(iii) The sides {± 7 tL} X (— 00 ,1) are mapped to the slit. 

The result now follows by taking the logarithm of the modulus of the map / 
in Equation (6.1). □ 

By using Proposition 6.1, we can obtain an explicit expression for the 
leading order wave velocity ci, and a Fourier series for the leading order 
surface profile 

Proposition 6.2 (ci and 77*). The leading-order wave velocity c\ and surface 
profile i]* are given by 

ci = coth(l/L), 

47 tL 

00 

r]* = ~a 2 T2 7 ~ 9 2 e~ n/L cos(nx/L), 

4tt 2 ' gL 2 + a 2 n 2 

n=l 

respectively. 


Proof. Recall how the wave velocity appeared on the right-hand side of 
Equation (3.3). By using the final part of Theorem A.l, we find 


ci = — 


f"(0) 


1 


coth(l/L), 


47r \ /'(0) ) 47 tL 

where / is the conformal map introduced in Equation (6.1) in the proof of 
Proposition 6.1. 

We now move to the surface profile. From [41] we know that 


( 6 . 2 ) ??* = -{g~a 2 dl) 1 

where x is defined by 


1 


CivL 


X 


2-jtL 


Xdp 


—ttL 


X(x) := ciG y (x, 1) + -G y (x, l) 2 . 


Written out, we have 


X(x) = 


1 cosh(l/L) cos {x/L) — 1 


87 r 2 L 2 (cos (x/L) — cosh(l/L )) 2 
with the elementary antiderivative 

#/ x = _ 1 _ sin(x/L) 

87 r 2 L cos(x/L) — cosh(l/L) ’ 
In particular, this means that 

/ 7TL 

Xdp = X # (ttA) - x B (-7tL) = 0, 

-7 tL 

so that Equation (6.2) reduces to 

(6.3) V* = ~(g - a 2 dinx- 
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In order to find the Fourier series for 77 *, we require the Fourier series of X- 
We may write 


X^x) = y 


1 + 


D —ix/L—l/L 


1 


^ g— ix/L—l/L ^ gix/L—l/L 

which, by expanding into geometric series, means that 


X \x) = i sgn{n)e-W/ L e inx / L . 


Hence, by termwise differentiation, we obtain 

^ OO 

X0) = 47 r 2 L 2 Y ne ~ n/L cos (nx/L), 

n =1 

which, combined with Equation (6.3), yields the result. □ 

One may note that 

ci = + 0(1/L 2 ) 

as L -> oo, which agrees with the speed of the solitary waves on infinite 
depth. When L is large, the surface profile is very similar to the surface in 
the localized case, see Figure 4b. At the other extreme, the first terms in the 
Fourier series will dominate. 


• 10“ 2 - 10” 2 



x 


(a) L = 1 (Whole period) 



x 

(b) L = 100 (Only part of a period) 


FIGURE 4. The leading order surface profile term, r/*, when 
g = l,a 2 = 0.01, cf. [41, Figure 1]. 


Appendix A. 

In this appendix, we provide two theorems that are used in order to get 
exact expressions for the rotational part of the stream function. Except for 
the final part, Theorem A.l is a standard result [35, p. 166]. Theorem A.2 is 
a less well known extension of Theorem A.l. 
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Theorem A.l (Green’s functions in M 2 ). Suppose that ft C M 2 is a simply 
connected domain and that zo £ ft. Furthermore, suppose that f: 12 —> D is 
a bijective conformal m.ap onto the open unit disk, extending continuously to 
a function ft —> D and satisfying f(zo) = 0. Then the function tp: ft 
defined by 

¥>(«) : = ^log(|/(z)|) 

is in extends continuously to the boundary of ft, and satisfies 

A</? = 5 Z0 , 

T\an = °- 


Furthermore, the harmonic function h defined by 


satisfies 


h(z) ■= ip{z) - ^ log (\z - zq\) 


X7U ) 1 (f"( Z 0 ) 

VW= *rlfS) 


after identifying M 2 and C via ( x , y) i— > x + iy. 


Proof. We first check the boundary values of the function ip. By assumption, 
/ extends continuously to dft, and every point on dft must necessarily be 
mapped to the unit circle. It is thus immediate that ip also extends continuosly 
to the boundary, and moreover, vanishes there. 

Identify now M 2 and C. Observe that since f(zo) = 0, we have 

f{z) = g(z)(z - z 0 ), zeft 

for some holomorphic function g, where \g\ > 0. Indeed, we must have 
g(zo) = f' ( zq ) ^ 0 because / is injective, and the injectivity of / also ensures 
that there can be no other roots. Thus 

ip(z) = ^ log (\z - z 0 |) + h(z), 

where 

h(z) := — Re log (g(z)) 

Z7T 

is harmonic by |^| > 0 and the Cauchy-Riemann equations. Hence, by 
Proposition 2.1, the function (p is L\ and satishes 

A ip = S Zo . 


The last assertion follows by observing that one must necessarily have 
g'(z 0 ) = lf"(zo), meaning that 


(^ lo g(9('))) 


(zo) 


i g'jzp) 

2vr g(z 0 ) 


1 f"(zo) 

4vr f'(zo) ’ 


whence we deduce from the Cauchy-Riemann equations that 


Vh(zo) 


1 (f"(z Q )\ 
4vr V f'(zo) ) 


□ 
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Theorem A.2 (Green’s functions in R 2 , mixed). Suppose that Q C R 2 is 
a simply connected domain and that zq £ Furthermore, assum.e that 
dfl = r n U where Tjy is C 1 and open in dfl. Finally, suppose that 
f: —> D \ ((—1, —a] x {0}), where a > 0, is a bijective conformal map of Q 
onto the unit disk with a slit, satisfying f(zo) = 0 and extending continuously 
to the boundary. This map should send Tp, to the unit circle and Tn to the 
interval (—l,a] x {0} ; and should extend analytically across Tjg (when viewed 
as a map on C). Then the function ip: ^ R defined by 

v( z ) ■= 7^ log (|/0)|) 

is in extends continuously to the boundary and satisfies 

Ap = S Z0 , 

P\r D = 

^ n T\r N = 11 ’ 

where d n denotes the normal derivative. 

Proof. The only change from Theorem A.l is checking that the normal 
derivative vanishes on r^r. This follows by using conformality. □ 
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